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ABSTRACT 

We present the detection of an intermediate-mass pre-main-sequence population embedded in the 
nebular filaments surrounding the 30 Doradus region in the Large Magellanic Cloud (LMC) using 
HST/NICMOS. In addition to four previously known luminous Class I infrared "protostars," the NIC- 
MOS data reveal 20 new sources with intrinsic infrared excess similar to Galactic pre-main sequence 
stars. Based on their infrared brightness, these objects can be identified as the LMC equivalent of 
Galactic pre-main sequence stars. The faintest LMC Young Stellar Objects in the sample have colors 
similar to T Tauri and have about the same brightness as T Tauri if placed at the distance of the LMC. 
We find no evidence for a lower-mass cut-off in the initial mass function. Instead, the whole spectrum 
of stellar masses from pre-main sequence stars with ~1.5M0 to massive O stars still embedded in dense 
knots appears to be present in the nebular filaments. The majority of the young stellar objects can be 
found to the north of the central starburst cluster R136. This region is very likely evolving into an OB 
association. 

The observations provide further evidence that star formation in the 30 Doradus region is very similar 
to Galactic star formation, and confirm the presence of sequential star formation in 30 Doradus, with 
present-day star formation taking place in the arc of molecular gas to the north and west of the starburst 
cluster. 

Subject headings: (galaxies:) Magellanic Clouds — galaxies: stellar content — stars: formation — 
stars: pre-main sequence 



1. INTRODUCTION 

The study of extragalactic star formation began 20 years 
ago with the identification of the first candidate protostar 
in N159 in the Large Magellanic Cloud (LMC) by Gatley 
et al. (1981). It was subsequently expanded by the iden- 
tification of other candidate "protostars" , in particular in 
the 30 Doradus region (Hyland et al. 1992). 30 Doradus 
is the most luminous giant HII region in the Local Group 
(Kennicutt 1984). The starburst region consists of multi- 
ple stellar generations ranging from 2 Myr old Wolf-Rayet 
and 03 main-sequence stars to evolved blue and red super- 
giants with ages up to 25 Myr (Walborn & Blades 1997; 
de Koter et al. 1998; Massey & Hunter 1998; Grebel & 
Chu 2000; and references therein). The stellar popula- 
tion has been studied with high spatial resolution with 
the Hubble Space Telescope (HST) in the optical down to 
2.8 Mq (Hunter et al. 1995) and with adaptive optics in 
the near-infrared (Brandl et al. 1996). Recently, Sirianni 
et al. (2000) extended the study of the initial mass function 
down to 1.35 M0. Infrared and radio observations reveal 



that molecular gas and warm dust is concentrated in an 
arc to the north and west of the central starburst cluster 
in 30 Doradus (Werner et al. 1978; Johansson et al. 1998; 
Rubio et al. 1998). Three early 0-type stars embedded 
in dense nebular knots (Walborn & Blades 1987, 1997), 
four luminous infrared "protostars" (Hyland et al. 1992) 
as well as numerous infrared sources associated with on- 
going star formation (Rubio et al. 1992, 1998) have been 
identified in the arc. Using HST/NICMOS (Thompson 
et al. 1998), we studied the embedded IR sources in the 
nebular filaments and their environment in a multi-color 
survey, aiming at a better understanding of the extent of 
triggered star formation, the nature of the embedded stel- 
lar population, and its luminosity function (see Walborn 
et al. 1999 for first results). The results of a complemen- 
tary HST/NICMOS program, which aimed at establishing 
a deep H-band luminosity function for stars in the central 
starburst cluster, are reported elsewhere (Zinnecker et al. 
1999, and in prep.). 

Section 2 gives an overview on the observations. Section 



^Based on observations with the NASA/ESA Hubble Space Telescope obtained at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy, Inc., under the NASA contract NAS5-26555, and observations at the European 
Southern Observatory, La Silla, obtained during technical time. 
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3 details the identification of pre-main sequence stars. In 
Section 4, the new results are discussed with respect to 
previous studies of young stellar objects in the LMC, and 
Section 5 presents a brief summary. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Ground-based near-infrared imaging 

Ground-based near-infrared imaging observations in the 
J-, H-, and K-band were carried out on 1993 February 15 
with the ESO/MPI 2.2m telescope and the IRAC-2a cam- 
era (Moorwood et al. 1992) at La Silla, Chile. IRAC-2a 
was equipped with a NICM0S3 array, and a pixel scale 
of 0'.'49 pixel"^ was chosen. The data were obtained as 
a 3x3 mosaic covering an area of approximately 350" x 
350". Integration times were 60s per position and filter, 
and the seeing was I'.'O to 1"2. 

The data were reduced in the standard way, used to 
identify embedded red objects, and — combined with the 
data obtained by Rubio et al. (1998) — helped to plan 
the HST/NICMOS observations. Figure |l| shows a 200" x 
225" subset of the K-band mosaic centered on 30 Doradus. 

2.2. Space-based near-infrared imaging 

Observations of 17 fields centered on infrared sources 
and one background field ("sky") were obtained with 
HST/NICMOS (NIC2) in Feb. and Mar. 1998. 14 of the 
fields were observed in the filters FllOW (320s), F160W 
(384s), and F205W (448s). Field 17 was only observed 
in selected narrow-band filters, i.e., F187N (Paa), F212N 
(H2), F215N (cont.), and F216N (Br7), while observations 
of fields 4 and 5 were only partially executed due to in- 
strument problems (NICMOS suspended due to a particle 
hit). The location and orientation of the individual NIC2 
fields are outlined in Figure |l|. In the present paper, we 
focus on the photometric analysis of the broad-band data. 

The data were reduced using synthetic darks customized 
for the appropriate detector temperature. The observa- 
tions of the background field in F205W were combined 
to produce a sky frame. This "sky" was subtracted from 
the on-source observation in F205W in order to remove 
the thermal emission of the telescope. Source identifica- 
tion and photometry were carried out with the IRAF im- 
plementation of DAOPHOT. The uncertainties quoted in 
Tables to | are DAOPHOT fitting uncertainties. Trans- 
formations from the HST system to ground-based photom- 
etry were computed by a comparison of HST/NICMOS 
observations of stars from the NICMOS list of stars for 
photometric transformations (NICMOS Data Handbook 
Version 4) to photometry of the same stars as presented 
in the 2MASS point source catalog (see also Skrutskie et 
al. 1997; Beichman et al. 1998). We derived the following 
color transformations (in units of [mag]): 

miio-J (0.263±0.020)*(niiio-mi6o) + (0.083±0.022) 

mi6o-H = (0.072±0.041)*(miio-mi6o)-l-(0.032±0.045) 

01205 -Ks = (0.268±0.063)*(mi6o-m205)-(0.038 ±0.070) 

The color transformations are derived for standard stars 
with J-H colors ranging from -0.2™ to -1-2.0" and H-Kg 
colors ranging from +0.1™ to +0.9™. The color transfor- 
mations from FllOW and F160W to ground-based J and 



H magnitudes are within the uncertainties in good agree- 
ment with the color transformations presented by Origlia 
& Leitherer (2000). Transformations from the HST NIC- 
MOS system to the CTI/CTIO photometric system, and 
possible complications due to spectral features which are 
not in common between standard stars and program stars, 
are discussed by Stephens et al. (2000). 

We estimate the absolute photometric uncertainty to be 
of the order of ±0™1. The brightness limit for sources de- 
tected in all three broad-band filters is mx ~ 19™8 for a 
signal-to-noise ratio of «7. Relative photometric uncer- 
tainties for individual sources are quoted in Tables |l| to 

i 

Positions of individual sources are based on the astro- 
metric information in the FITS file headers. The rela- 
tive astrometric uncertainty in the position of individual 
sources within one NICMOS frame should be of the or- 
der of ±0"02, whereas the absolute positional uncertainty 
might be as large as ±1". 

3. IDENTIFICATION OF PRE-MAIN SEQUENCE STARS 

Near infrared (NIR) color-color diagrams can be used 
to separate stars with infrared excess due to extinction 
from stars with intrinsic infrared excess. Intrinsic infrared 
excess can be attributed to the presence of circumstellar 
material, and thus is a sign of youth. As pointed out by 
Lada & Adams (1992), different types of objects occupy 
well-defined, distinct regions in NIR color-color diagrams. 
Many young stellar objects lie outside the normal redden- 
ing band. Infrared "protostars" (Class I sources) typically 
are highly reddened (Ay = 10™ to 40™) and have both J~H 
and H-K colors > +1™5. Herbig AeBe stars tend to have 
J-H and H-K colors between +0™8 and +1™5. Classical 
(CTTS) and weak-line T Tauri stars (WTTS) occupy a re- 
gion characterized by lower extinction (Ay < 10™). CTTS 
cluster around J-H « +1™0 and H-K w +0™8, whereas 
WTTS typically fall within the reddening band. Classical 
Be stars on the other hand can be found near the blue end 
of the main-sequence. They show less intrinsic reddening 
than Herbig AeBe or T Tauri stars and tend to have J-H 
colors between — 0™2 and +0'."3 and H-K colors between 
-0™4 and +0'.°3. 

Variability of individual sources and other effects like, 
e.g., source geometry make it difficult to determine the 
physical properties of individual young stellar objects. The 
evolutionary nature of a population of young stellar ob- 
jects can, however, be deduced on a statistical basis from 
their location in a J-H vs. H-K color-color diagram. Ta- 
bles and H list all sources which were detected in all three 
bands and which exhibit intrinsic IR excess (i.e., they lie 
within the photometric uncertainties clearly outside the 
standard reddening band). 

Figure ^, top, shows NIR color-color diagrams of two 
fields centered approximately 30" and 50" to the west of 
R136a. In addition to stars with intrinsic colors close to 
the colors of unreddened dwarf and giant stars, some stars 
show clear evidence of IR excess. The majority of the 
stars with IR excess are located in the area of the normal 
reddening band, and their IR excess is very likely due to 
enhanced (foreground) extinction of Ay = 5™ to 10™. 

Figure g, bottom, shows the color-color diagrams for two 
of the fields located to the north of R136a. These fields 
(30Dor-NIC 11, 12, 13, 14, and 15) house a relatively large 
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Fig. 1.— 200" X 225" K-band image (obtained with IRAC-2a at the MPI/ESO 2.2m telescope) of the 30 Doradus region indicating the 
location and orientation of the HST/NICMOS (NIC2) pointings. North is up and East is to the left. 



Table 1 

Candidate Class I sources and Herbig AeBe/TTauri stars to the south and west of R136a 



Name 


Alias 


a(2000) 


(5(2000) 


J 


J-H 


H-K, 








[hms] 


[ 


O f 


"] 


[mag] 


[mag] 


[mag] 


30Dor-NIC08a 




5 


38 29.87 


-69 


06 


11.3 


16.95±0.06 


0.44±0.07 


0.55±0.14 


30Dor-NIC07a 




5 


38 33.14 


-69 


06 


11.1 


18.68±0.12 


1.84±0.15 


1.75±0.11 


30Dor-NIC07b 




5 


38 33.78 


-69 


06 


15.4 


20.03±0.11 


0.86±0.14 


0.79±0.14 


30Dor-NIC03a 


[HJ91] P3\ IRSW-302 


5 


38 34.57 


-69 


05 


58.2 


15.23±0.09 


2.16±0.10 


2.05±0.05 


30Dor-NIC02a 


5 


38 36.72 


-69 


06 


21.0 


18.40±0.09 


0.83±0.11 


0.78±0.08 


30Dor-NIC01a 




5 


38 40.12 


-69 


06 


37.4 


19.28±0.09 


0.36±0.13 


0.65±0.17 


30Dor-NIC01b 




5 


38 40.21 


-69 


06 


33.2 


19.21±0.07 


0.78±0.09 


0.67±0.06 


30Dor-NIC09a 




5 


38 43.47 


-69 


07 


00.7 


19.97±0.10 


0.62±0.13 


0.58±0.13 



References.— 1 Hyland & Jones (1991), ^ Rubio et al. (1998) 
number of stars with intrinsic NIR excess. Based on their 



locations in the NIR color-color diagram, the objects ap- 
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Table 2 

Candidate Class I sources and Herbig AeBe/TTauri stars to the north and north-east of R136a 



Name 




Alias 






a(2000) 

[hms] 


(5(2000) 
[° ' "] 






J 

[mag] 






J-H 
[mag] 






H-Ks 
[mag] 


30Dor- 


-NIC16a 


[HJ91] 


P2^ 




5 


38 41.62 


-69 


03 


54, 


.7 


15, 


,38±0, 


,05 


2, 


.11±0, 


,10 


1, 


,56±0. 


,12 


aODor- 


-NICiaa 






5 


38 44.65 


-69 


05 


01, 


.6 


18, 


,03±0, 


,06 


0, 


.60±0, 


,07 


0, 


,57±0. 


,05 


aODor- 


-NlCiab 








5 


38 45.00 


-69 


04 


57, 


.9 


18, 


,90±0, 


,06 


0, 


.71±0, 


,08 


1, 


,13±0. 


06 


aODor- 


-Niciac 








5 


38 45.10 


-69 


04 


54, 


.9 


19, 


,46±0, 


,08 


0, 


.74±0, 


,11 


0, 


,79±0. 


,09 


aODor- 


-NIClla 








5 


38 45.96 


-69 


05 


17, 


.0 


20, 


,05±0, 


,13 


0, 


,62±0, 


,17 


0, 


,74±0. 


,14 


aODor- 


-NICllb 








5 


38 46.57 


-69 


05 


20, 


,1 


17, 


,65±0, 


,04 


0, 


,64±0, 


,05 


1, 


.22±0. 


,04 


aODor- 


-NIC12a 








5 


38 46.99 


-69 


05 


08, 


.3 


19, 


,46±0, 


,08 


1, 


.05±0, 


,10 


1, 


,55±0. 


,07 


aODor- 


-NIC12b 


[HJ91] 


Pl\ 


IRSN-1222 


5 


38 47.11 


-69 


05 


06, 


,1 


17, 


,46±0, 


,09 


2, 


,16±0, 


,11 


2, 


,09±0. 


,08 


aODor- 


-NIC12C 






5 


38 47.24 


-69 


05 


03, 


.4 


19, 


,86±0, 


,12 


1, 


.34±0, 


,14 


1, 


,09±0. 


,12 


aODor- 


-NIC12d 


[HJ91] 


Pl\ 


IRSN-1262 


5 


38 47.25 


-69 


05 


02, 


.3 


17, 


,15±0, 


,06 


2, 


.83±0, 


,07 


2, 


,54±0. 


,04 


aODor- 


-NIC15a 






5 


38 48.30 


-69 


04 


10, 


.3 


19, 


.01±0, 


,08 


0, 


.77±0, 


,10 


0, 


,67±0. 


,09 


aODor- 


-NIC15b 


[HJ91] 


P4\ 


IRSN-I342 


5 


38 48.45 


-69 


04 


12, 


.0 


17, 


,06±0, 


,07 


1, 


,82±0, 


,08 


1, 


,76±0. 


,04 


aODor- 


-NIC14a 






5 


38 48.76 


-69 


04 


47, 


.6 


19, 


,16±0, 


,09 


0, 


,61±0, 


,11 


0, 


,64±0. 


,12 


aODor- 


-NIC14b 








5 


38 48.82 


-69 


04 


47, 


,1 


19, 


,91±0, 


,12 





,92±0, 


,16 


0, 


,97±0. 


,17 


aODor- 


-NIC12e 








5 


38 48.87 


-69 


05 


02, 


.2 


20, 


,30±0, 


,14 





.79±0, 


,18 


0, 


,89±0. 


,16 


aODor- 


-NIClOa 








5 


38 55.55 


-69 


05 


25, 


.8 


16, 


.35±0, 


,03 


0, 


.36±0, 


,04 


0, 


,41±0. 


03 



References.— 1 Hyland & Jones (1991), ^ Rubio et al. (1998) 



pear to be Class I "protostars" and Herbig AeBe stars or 
TTauri stars (see Lada & Adams 1992). Figure |^ shows 
a color-magnitude diagram for the objects listed in Tables 
1 to 3. The faintest stars in the sample have 3Ri2{F^0. 
According to the 2MASS point source catalog (Second In- 
cremental Data Release, Cutri et al. 2000), TTauri itself 
has J=7^26, J~H=1™02, and H-K=0"^91. If moved from 
its location in the Taurus T association (distance 140 pc) 
to the distance of the LMC (distance w 50 kpc), TTauri 
would have J=20™0. Thus the faintest objects listed in Ta- 
ble and ^ could actually be LMC counterparts to T Tauri. 

The other fields to the west of the cluster (30Dor-NIC 
7, 8) have very few stars with intrinsic NIR excess. The 
same is the case around 30Dor-NIC 10, which is located 
NE of the cluster within the arc, and NIC 6 and 9, which 
are located to the south and south-west of the cluster. 

Thirteen of the fourteen fields with broadband obser- 
vations in all three bands house at least one source with 
intrinsic NIR excess (see Tables |l| and ^) . Twelve of these 
fields are shown in Figure |^ with the NIR excess sources 
identified. 

The number of candidate sources identified is very likely 
a lower limit to the total number of young stars in the 
nebular arc as we required candidates to be detected in 
all three bands, and as not all young stellar objects ex- 
hibit strong intrinsic IR excess. With six detections, the 
30DOR-NIC12 field exhibits the highest incidence of in- 
trinsic IR excess sources. This is probably due to the fact 
that the location of 30DOR-NIC12 coincides with one of 
the density peaks ("pillar") in the arc. 

The comparison with ground-based photometry is re- 
stricted to the objects which are in common between our 
study and the studies by Hyland et al. (1992) and Rubio et 
al. (1998): the four "protostars" PI to P4 and four objects 
resolved in knots 1 to 3. The aperture size of 5" as used 
by Hyland compared to our diffraction limited resolution 



of O'.'l (FllOW) to 0'.'2 (F205W) makes a comparison not 
very meaningful. The seeing limited observations by Ru- 
bio et al. (1998) have a better resolution of wl", but are 
still by about a factor of 6 to 10 lower spatial resolution 
than the HST NICMOS observations. Hence a number 
of fainter background sources contribute to the integrated 
photometry presented by Rubio et al. (1998), which are 
not included in our photometry due to the higher spatial 
resolution. 

A direct comparison to the photometry presented by 
Rubio et al. (1998) does not yield any evidence for large, 
systematic differences. The scatter in the photometric 
measurements for individual sources can be explained by 
the differences in spatial resolution, the non-standard HST 
NICMOS passbands (in particular F205W), and last but 
not least by intrinsic variability of individual YSOs. 

4. COMPARISON TO PREVIOUS WORK ON YOUNG 
STELLAR OBJECTS IN THE 30 DORADUS REGION 

4.1. Early 0-type stars embedded in nebular knots 

Based on optical spectroscopy, Walborn & Blades (1987, 
1997) identified three early 0-type stars associated with 
nebular knots (knots 1 to 3, see Figure ^ in the north- 
east and the west of 30 Doradus as a first evidence for 
ongoing star formation in the 30 Doradus region. Similar 
to knot 2, which was originally resolved into individual IR 
sources by Rubio et al. (1992), the HST/NICMOS obser- 
vations revealed also knots 1 and 3 as compact multiple 
systems, very much like the Trapezium system in the Orion 
cluster (Walborn et al. 1999). Photometry and positions of 
the brightest components of the multiple systems are sum- 
marized in Table 0. The near infrared colors corroborate 
the identification of the subcomponents as early-type stars 
seen through a moderate amount of foreground extinction. 

Assuming that the ratios of projected separations are 
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H - K [mag] H - K [mag] 



Fig. 2. — Examples of necir-infrared color-color diagrams of NIC2 fields to the west (top) and to the north of R136a (bottom). The location 
of (unreddened) main-sequence stars is indicated by a solid line, the location of giant stars by a dashed line. The dotted lines indicate the 
reddening band for normal extinction. Asterisks mark steps of 5™ in visual extinction. Stars with intrinsic infrared excess fall outside the 
region of normal reddening. In addition to previously known luminous Class I sources, stars with the color characteristics of intermediate- to 
low-mass young stellar objects (Herbig AeBe and TTauri stars) are detected. 



Table 3 

Near infrared photometry of the brightest individual objects in knots 1 to 3 



Name 


Alias^ 


a(2000) 
[hms] 


5(2000) 
[° ' "] 


J 

[mag] 


J-H 
[mag] 


H-Ks 
[mag] 


Knot la 


IRSN-101 


5 38 45.37 


-69 05 08.9 


14.67±0.08 


0.19±0.09 


0.03±0.05 


Knot lb 


IRSN-96 


5 38 45.20 


-69 05 08.2 


15.34±0.06 


0.22±0.07 


0.20±0.06 


Knot Ic 


IRSN-96 


5 38 45.11 


-69 05 08.1 


16.63±0.10 


0.38±0.14 


0.58±0.19 


Knot 2a 


IRSN-135 


5 38 48.05 


-69 04 43.5 


14.82±0.06 


0.31±0.07 


0.07±0.04 


Knot 2b 




5 38 48.00 


-69 04 43.5 


17.01±0.09 


0.31±0.12 


-0.14±0.14 


Knot 2c 




5 38 47.96 


-69 04 43.6 






2 


Knot 3ab^ 


IRSW-29 


5 38 34.55 


-69 06 07.2 


15.61±0.17 


0.23±0.24 


0.01±0.21 


Knot 3c 




5 38 34.51 


-69 06 06.9 


16.99±0.05 


1.15±0.07 


0.28±0.05 



Note. — ^ Rubio et al. (1998), ^ mx = 16™28±(K"08, * close binary, photometry of unresolved system, coordinates for center of light 



representative for their true separations, we find the fol- 
lowing ratios (separation of close pair divided by separa- 
tion to third component): 

• Knot 1: 1:2.7 (O.lSpc vs. 0.35 pc) 

• Knot 2: 1:1.8 (0.054 pc vs. 0.097pc) 

• Knot 3: 1:2.8 (0.033 pc vs. 0.090 pc) 

For comparison, the four central stars in the Trapez- 
ium cluster in Orion have projected separations between 



«0.02 pc and w0.04 pc. Thus similar to the Trapezium sys- 
tem, the three early-type triple systems in the 30 Doradus 
Nebula appear non-hierarchical. Since non-hierarchical 
multiple systems are dynamically unstable (e.g, Mirzoian 
& Salukvadze 1985; Stcrzik & Durisen 1998; Orlov & 
Petrova 2000; and references therein), the close associa- 
tions of early-type stars in the 30 Doradus Nebula provides 
further evidence for their youth. 
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Fig. 3. — J vs. J-Ks color-magnitude diagram for candidate Young Stellar Objects and the brightest objects in Knot 1 to 3 (see Tables 1 
to 3). The approximate J-band brightness for subsamples corresponding to candidate Class I sources, Herbig AeBe stars and T Tauri stars 
in the LMC is indicated. 



As already noted by Walborn et al. (1999), the multiple 
early-type systems in knots 1 and 2 are associated with 
a larger number of fainter stars, and might thus consti- 
tute very young open clusters. A third cluster of faint, 
extincted stars can be found around 30Dor-NIC15b (see 
Figure I and I). The on average somewhat larger physical 
separation of the early type stars in the center of these 
open clusters compared to the Orion cluster might indi- 
cate that the former already started to disperse gradually 
(see, e.g., Kroupa et al. 2001). Alternatively, the physical 
conditions which led to the formation of the 30 Doradus 
young open clusters might have been different (less ex- 
treme?) than the conditions which led to the birth of the 
Orion cluster. A more detailed discussion of the stellar 
populations in the association and the young open clus- 
ter candidates will be presented in a forthcoming paper 
(Grebel et al., in prep.). 

4.2. Evidence for a low-mass cut-off among pre-main 
sequence stars in 30 Doradus? 

In their near infrared survey of the 30 Doradus region, 
Hyland et al. (1992) found four infrared candidate "proto- 
stars" with K magnitudes brighter than 12™3, and masses 
in the range of 15 to 20 Mq (as deduced from the 
bolometric luminosities of the infrared sources). Based 
on the lack of objects with strong infrared excess and 
K-magnitudes between 12™ 3 and 13™ 3 (the latter being 
the brightness limit of their survey), Hyland et al. (1992) 



suggested that there might be a lower-mass cut-off below 
15 M0 in the initial mass function for the present-day star 
formation in 30 Doradus. 

The color-magnitude diagram (Figure ^) reveals that 
the population of sources with intrinsic IR excess extends 
down to the sensitivity limits of our observations. The 
candidate pre-main sequence stars listed in Tables |] & 
Hhave near-infrared colors (and considering the distance 
to the LMC, brightness) quite similar to Galactic Her- 
big AeBe and T Tauri stars with masses in the range from 
«1.5M0 to «7M0. Furthermore, [HJ 91] PI is resolved 
into two red objects (30Dor-NIC12b and d, see Figure |^) 
with the fainter one having mx = 13™2 (see also Rubio et 
al. 1998). Thus our HST/NICMOS survey gives no evi- 
dence for a cut-off in the present-day initial mass function 
possibly down to at least 1.5 Mq assuming the infrared ex- 
cess sources have properties similar to Galactic pre-main 
sequence stars. 

4.3. The low-mass IMF in the central starburst cluster 

Sirianni et al. (2000) analyzed archived WFPC2 data of 
the 30 Dor region. The identification of pre-main sequence 
stars with masses down to 1.35 M© close to the center of 
the R136 cluster by Sirianni et al. (2000) is based on the 
assumption that there is very little differential reddening 
towards 30 Doradus. Only then can pre-main sequence 
star candidates be identified exclusively by their position 
in a color magnitude diagram. 
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Fig. 4. — Finding charts for prc-main sequence candidates in 12 of the NICMOS/NIC2 fields in 30 Doradus. The field of view of each 
individual field is 19" x 19", and the images shown have been observed through the F205W filter. Note that 30Dor-NIC12f is identical to 
30Dor-NIC13a. The locations of knot 1 to 3 (Walborn & Blades 1987, 1997) are indicated by circles. 



Of the NIC2 fields studied by us, only the fields to the 
west of R136a overlap with the HST wide field CCD data 
(though not the planetary camera data) studied by Siri- 
anni et al. (2000). The NICMOS data clearly reveal the 
presence of substantial differential redding in the overlap- 
ping fields. Variable extinction close to the cluster cen- 
ter was reported by Brandl et al. (1996) in their multi- 



wavelength study of the initial mass function of the star- 
burst cluster, and is also evident in the data presented by 
Rubio ct al. (1998) and Scowcn et al. (1998). This indi- 
cates that a certain fraction of the purportedly prc-main- 
sequcncc stars in the HST/WFPC2 sample by Sirianni et 
al. (2000) might actually be main-sequence stars viewed 
through 5™ to 10™ of visual extinction. An accurate ini- 
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tial mass function can hence only be derived by applying 
extinction corrections to individual sources, which in turn 
requires multi-color photometry or spectroscopy of indi- 
vidual stars (see also, e.g., Selman et al. 1999; Panagia et 
al. 2000). 

The scarcity of sources with intrinsic infrared excess in- 
side the 30 Dor nebular arc, however, does not rule out 
the presence of more pre-main sequence stars. The early- 
type stars in the 30 Doradus cluster create an intense ra- 
diation field, which has an adverse effect on circumstcUar 
disks. Observations of Galactic HIT regions like Orion or 
NGC 3603 revealed that externally illuminated circumstel- 
lar disks get photo-evaporated and disperse on short time 
scales of typically 10^ to 10^ yr (e.g., Henney & O'Dell 
1999, Brandner et al. 2000a, RichUng & Yorke 2000). 

4.4. Overluminous Herbig AeBe stars? 

Based on data obtained by the EROS project, Lamers 
et al. (1999) identified a group of irregular variable stars in 
the bar of the LMC with spectral characteristics similar to 
Galactic Herbig AeBe stars. The majority of these stars, 
however, appear to be more luminous than Galactic Her- 
big AeBe stars, and are located above the birthline (Palla 
& Stabler 1993) of Galactic Herbig AeBe stars. Lamers 
et al. (1999) argue that the higher luminosity might be 
either due to higher overall accretion rates or due to the 
fact that pre-main sequence stars in the LMC become opti- 
cally visible at an earlier phase than their Galactic equiva- 
lents (smaller dust-to-gas ratio in the LMC). Alternatively, 
Lamers et al. (1999) propose that they might detect the 
high end of the luminosity function of pre-main sequence 
stars in the LMC. 

In order to test if the infrared luminous "protostars" 
from the work by Hyland et al. (1992) could be rep- 
resentatives of this group of "overluminous" pre-main- 
sequence stars, we compare their near infrared photometry 
to the Galactic Herbig AeBe star Z CMa. Z CMa has a J- 
magnitude of 5™79, a J-K color of 2™25, and is located at 
a distance of «1150pc (e.g., Leinert et al. 1997). If placed 
at the distance of the LMC, ZCMa would have mj=14™0 
and mK=ll™7, i.e. quite in the range of brightness values 
observed for the "protostars" from Hyland et al. (1992). 
The redder J-K colors of the candidate "protostars" iden- 
tified by Hyland et al. (1992) compared to ZCMa could 
be due to differences in the geometry of circumstellar ma- 
terial. Edge-on circumstellar disk sources in Galactic star 
forming regions can have J-K colors ^7™ (e.g. Brandner et 
al. 2000b, Table 1). Thus again we find no strong evidence 
that star formation in 30 Doradus might be different from 
Galactic star formation. This is also in agreement with the 
study of pre-main sequence candidates around SN 1987A 
in 30 Dor C by Panagia et al. 2000, who do not find ev- 
idence for a large population of overluminous pre-main 
sequence stars (e.g., Fig. 4 in their paper). 



These findings do not exclude, however, that star forma- 
tion in lower mctallicity environments, like the bar region 
of the LMC, might be different. 

5. SUMMARY 

We analyzed HST/NICMOS data and groundbased in- 
frared imaging obtained with the ESO/MPI 2.2m tele- 
scope and IRAC2 of the 30 Doradus region. We find clear 
evidence for ongoing, extended star formation in a 10 pc x 
15 pc region to the north of the central starburst cluster. 
This region is part of the arc of molecular gas and warm 
dust around 30 Doradus. Star formation in this region 
may have been triggered by the starburst cluster. The 
fainter IR-excess sources have colors and brightness simi- 
lar to Herbig AeBe and T Tauri stars with masses possibly 
in the range of 1.5 M© to 7M0, and could thus be LMC 
equivalent of Galactic low- to intermediate-mass pre-main 
sequence stars. The present-day initial mass-function ap- 
parently shows no cut-off down to the mass (resp. mag- 
nitude) limit of our study. The whole spectrum of stellar 
masses from pre-main sequence stars with masses ^1.5 Mq 
to massive O stars still embedded in dense knots might be 
present in the nebular filaments. 

The spatial extent of pre-main-sequence stars in the re- 
gion to the north of 30 Doradus suggests that we are wit- 
nessing the birth of an OB association including a small 
number of open clusters. Similar older associations, in 
part with low-density clusters, are seen throughout the 30 
Doradus region (e.g., Walborn & Blades 1997; Grebel & 
Chu 2000). Star formation to the west and south appears 
to be less intense than in the large scale, nascent OB asso- 
ciation to the north of the cluster. This might be indica- 
tive of a different, more isolated mode of star formation. 
Present day star formation activity in 30 Doradus coin- 
cides with the spatial distribution of the densest features 
of the molecular gas, or their interfaces with the central 
cavity being evacuated by the R136 cluster. 

Follow-up near infrared spectroscopy of the pre-main 
sequence candidates is required in order to confirm their 
youth and to study their physical properties, such as spec- 
tral types or accretion rates. 
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